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1.  INTRODUCTION 


In  this  article  we  shall  deal  with  the  Wilson  Cloud  formed  by  low 
altitude  nuclear  explosions^.  Despite  the  fact  that  it  is  a well-known 
phenomenon,  we  have  found  no  detailed  treatment  of  it  in  the  literature. 

The  basis  of  the  effect  is  straightforward.  During  the  weak  shock  phase 
of  an  explosion,  over  a certain  volume,  the  pressure  behind  the  shock  front 
becomes  less  than  ambient  atmospheric  pressure.  Thus  parcels  of  air  passing 
through  the  shock  undergo  a sudden  expansion  and  on  occasion  temporarily 
drive  the  relative  humidity  to  supersaturation.  Under  these  conditions  the 
radii  of  water  droplet  aerosols  present  can  grow  many  orders  of  magnitude, 
thereby  momentarily  forming  a cloud  in  a region  behind  the  advancing  shock. 
We  shall  restrict  the  discussion  here  to  maritime  aerosols  composed  of  salt 
solution  droplets. 

We  shall  demonstrate  that  (1)  there  is  a threshold  of  ambient  rela- 
tive humidity  for  significant  cloud  formation  independent  of  the  explosion 
energy;  and  (2)  over  a wide  range  of  explosion  energies  the  maximum  cloud 
droplet  sizes  are  nearly  independent  of  the  ambient  aerosol  spectrum  and 
may  be  scaled  with  the  explosion  energy. 

In  Section  2,  below,  we  discuss  the  hydrodynamic  flow  behind  the 
shock  wave  and  its  effect  on  the  relative  humidity  which  drives  the  cloud 
formation.  In  Section  3 we  review  droplet  growth  kinetics  and  the  basis 
for  droplet  size  scaling.  Finally,  we  give  in  Section  4,  a universal 
scaled  formulation  of  droplet  size  and  cloud  location  together  with  a dis- 
cussion of  the  limitations  on  our  approach. 
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2.  SHOCK  HUMIDITY  PROFILES 


We  shall  concern  ourselves  here  with  the  time,  t,  history  of  rela- 
tive humidity,  S,  and  temperature,  T (°K),  of  an  air  parcel  passing  through 
the  shock  wave.  The  parcel  shall  be  denoted  by  its  initial  altitude,  z, 
and  its  radial  distance  r,  from  the  explosion  center.  S is  given  by 

S = [X/(X  + e)]  P/ps  (2.1) 

where  P is  the  pressure;  X is  the  absolute  humidity  (grams  of  water  vapor 
per  gram  of  dry  air);  e is  the  specific  gravity  of  water  vapor  relative  to 
that  of  air  (0.621);  and  P$(T)  is  the  saturation  vapor  pressure 

p (T)  = (3.53  x 104  dynes  cm'2)  exp  [-(Le/Rd)(l/T  - 1/300)]  (2.2) 

L is  the  heat  of  vaporization  of  water  (2.42  x lO1^1  erg  g’^);  RH  is  the 

-1  -1  ^ 
dry  air  gas  constant  (2.87  x 10  erg  g'  K ). 

A crucial  simplifying  approximation  results  from  neglecting  the 

water  condensation,  AX,  compared  to  X.  X is  typically  of  the  order  of 
_2 

10  and  for  explosion  energies  of  the  order  of  10  kt>AX  is  at  worst  1% 
of  X.  (We  assume  throughout  a typical  maritime  aerosol  number  density 

_3 

of  100  cm  ).  The  heat  of  condensation  associated  with  the  AX  gives  a 
temperature  deviation  from  adiabaticity 

AT  » -(L/C  )AX  (2.3) 

where  Cp  is  the  specific  heat  of  air  at  constant  pressure  (1.05  x 107 
erg  g"1  K'1).  Using  AX-10'4  we  find  AT-0.2K  and  (AS/S)  = -17.45  (AT/T) 

+ (AX/X)  ~1%  which  represents  the  strength  of  the  coupling  between  dry 
air  shock  hydrodynamics  and  cloud  formation.  Neglecting  this  coupling 
allows  the  use  of  existing  dry  air  calculations  of  shock  wave  profiles 
for  pressure  and  temperature.  These  calculations  assume  that  the  post 
shock  wave  flow  is  adiabatic,  i.e., 

T = Ts(P/Ps)(y'1)/y;  y = 1.4  (2.4) 
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where  P<.  and  are  the  pressure  and  temperature  of  the  parcel  immediately 
following  the  passage  of  the  shock  front.  Energy  scaling  laws  presented 
below  indicate  the  limitations  on  the  decoupling  approximation. 

Decoupling  of  the  dr/  air  shock  flow  from  the  cloud  formation  allows 
us  to  write 

S (r,  z,  t)  = SQ(z)  M(r,  z,  t)  (2.5) 


where  p 

M = (P/PQ)  exp  {17.45  -j-  -D)  (2-6) 

PQ>  T , PQ»  SQ  are  respectively,  the  ambient  pressure,  temperature,  density, 
and  relative  humidity  of  the  parcel  in  question.  M depends  on  (P/P  ) and 
(p/pQ).  These  ratios  are  functions  of  £ = (r/rQ)  and  t = (t/t  ) only. 

Here  rQ  and  t are  the  dynamic  length  and  dynamic  time  of  the  explosion^ : 


= (E/PJ1/3 


and 

t i E1/3P  '5/6  p 1/2 
O 0 Mo 

E is  the  hydrodynamic  energy  of  the  explosion.  Figure  1 gives  plots  of 

the  universal  function  M (t,  a)  versus  x for  various  £ values  (T  = 288.15 K 

is  used).  These  curves  are  based  on  the  numerical  treatments  of  shock 

(3) 

profiles  due  to  Needham,  et  al'  ' . Their  work  is  in  apparent  agreement 
with  that  reported  in  Reference  2 and  Reference  4. 

As  we  shall  see  below  significant  cloud  formation,  i.e.,  large  drop- 
let growth,  results  only  when  the  relative  humidity  in  a parcel  at  £ is 
driven  to  supersaturation  S > 1.0.  Thus,  from  Equation  (2.5)  there  is 

an  energy  independent  threshold  of  ambient  humidity  S^(£)  for  cloud 
formation  at  each  £ given  by  M*(£)’^  = (max{M(x ,£) })_1  . Figure  2 plots 
Sth(£)  versus  £.  (Ax(£)  is  defined  after  Equation  4.2).  We  find  from 
Figure  1,  that,  regardless  of  S0,  cloud  formation  can  begin  only  after 
time  t = 0.35  and  outside  a radius  £ - O.b.  We  emphasize  here 
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that  cloud  foimation  first  begins  during  a period  when  neither  strong  shock 
(self-similar)  theory  nor  weak  shock  (asymptotic)  theory  is  valid.  From 
x ~ .03  to  t ~ 1 a numerical  treatment  of  shock  hydrodynamics  is  required. 
We  note  in  passing  that  in  the  region  of  cloud  formation,  spatial  motion 
of  the  parcel  due  to  shock  winds  is  generally  small  and  of  no  interest  to 
the  problem.  A striking  feature  apparent  from  Figure  2 is  that  if  the 
ambient  humidity  is  nowhere  greater  then  70%  then  no  Wilson  cloud  of 
significance  can  form. 
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The  equation  of  growth  for  a salt  solution  droplet,  the  principal 

(51 

constituent  of  maritime  aerosols,  is  ' 


r(dr/dt)  = [S(t)  - Seq(r,m)]/Y(T) 


(3.1) 


where  S is  the  shock  driven  relative  humidity  in  the  air  parcel  contain- 
ing the  droplet  as  discussed  above.  Seq  (r,m)  is  the  equilibrium  relative 
humidity  for  a droplet  with  radius  r,  and  salt  mass  m.  Y(T)  is  a weakly 
temperature-dependent  coefficient  determined  by  the  heat  and  water  vapor 
diffusion  rates  in  air.  In  c.g.s.  units 

Seq(n,m)  = exp  (2cre/p^RdTr)  - 3i  mM/np^Wr3  (3.2) 

Y(T)  = (LpL/KT)(Le/(RdT)  - 1)  + PLRdT/(DePs(T) ) 

(Y  (300K)  = 7.26  x 105  sec/cm2)  (3.3) 

where  quantities  not  given  previously  are  defined  here: 

_3 

p,  = density  of  liquid  water  = 1 g cm 

L _i 

a = surface  of  water  = 72  dyne  cm 

3 -1-1 

K = thermal  conductivity  of  air  = 2.5  x 10  erg  cm  sec 

-1  2 -1 

D = water  vapor  diffusivity  in  air  = 2 x 10  cm  sec 

i = Van  Hoff's  constant  =2.7 

M = molecular  weight  of  water  = 18 

W = molecular  weight  of  salt  = 58 

As  an  illustration  of  the  droplet  growth  and  evaporation  cycle. 

Figure  3 considers  five  exemplary  droplet  species  of  salt  mass  mi  (grams) 
and  relative  population  (as  suggested  by  Reference  5)  contained  in  a 
typical  shocked  air  parcel.  (For  this  case  E is  of  the  order  of  10  kt  and 
SQ  ~ 0.89).  The  relative  humidity  S seen  by  the  dorplets  is  plotted  ver- 
sus their  radii  starting  from  their  entrance  to  the  shock  until  their 
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evaporation  behind  the  advancing  clouds.  The  S (r,  ) curves  are 

also  shown.  Upon  entering  the  shock  wave  a droplet  quickly  evaporates 

and,  after  a sudden  drop  in  S,  its  nucleus  returns  to  its  equilibrium 

curve.  It  remains  there  until  S is  driven  beyond  supersaturation,  at 

which  point  rapid  droplet  growth  begins.  Each  species  tagged  by  its 

initial  mass  is  found  to  grow  nearly  to  the  same  maximum  radius.  Such 

narrowing  of  the  droplet  radial  spectrum  is  in  common  with  droplet  growth 
f5l 

m thunderheads'  ' . This  feature,  which  we  shall  call  universal  growth, 
is  particularly  striking  in  Wilson  cloud  growth  where  the  droplet  size  at 
maximum  growth  are  independent  of  the  ambient  aerosol  spectrum. 

Universal  growth  may  be  understood  by  noting  that  for  each  droplet 
species,  once  S passes  above  unity,  where  significant  growth  results, 

S (r,m)  for  the  resulting  r is  very  nearly  equal  to  unity  and  may  be  so 
replaced  in  Equation  3.1.  This  allows  the  growth  equation  to  be  rewritten 
in  terms  of  scaled  variables  and  the  universal  function  M(x,£)  driving  the 
relative  humidity.  We  have  then,  independent  of  species,  i.e.,  independent 
of  m. 


where 


d( ( 1/2)R ^(x,£)/dx  = (SQ(z)  M(x,£)  - 1) 


r2(t,r,m)  = R2(t,£)  tj-1 


(3.4) 


(3.5) 


Y is  nearly  constant  over  the  course  of  growth  and  evaporation  and  we  use 
Y(300)  as  given  above.  At  a given  £, cloud  formation  begins  and  maximum 
growth  is  obtained  at  the  times  for  which  M(x,£)  = S'1.  We  denote  these 
times  as  x^*(£)  and  t2*U)  respectively.  They  may  be  easily  read  from 
Figure  1.  Since  droplets  typically  attain  radial  growths  significantly 
larger  than  their  ambient  radii,  a suitable  initial  condition  for  Equation 
3.4  is 


R2  (t1*(£),£)  = 0 


(3.6) 


independent  of  species.  No  new  length  scale  is  introduced  thereby  and  the 
solution  to  Equations  3.4  and  3.6  is  therefore  universal. 

To  clarify  the  meaning  of  universal  scaled  droplet  growth  (Equations 
3.4,  3.5  and  3.6)  note  that  they  imply  the  following:  if  at  a distance  r 
and  time  t from  an  explosion  of  energy  E the  droplet  radius  in  the  Wilson 
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cloud  is  r,  then  for  explosion  energy  E',r'  = r ( E ' /E) 1 at  r1  = r(EVE)1^3 
and  t'  = t (E1  / eP  assuming  the  same  ambient  humidity  prevails.  Further- 
more the  approximate  cloud  droplet  size  can  be  calculated  entirely  without 
regard  to  the  ambient  aerosol  spectrum. 


T» 


4.  RESULTS  OF  SCALING  AND  UNIVERSAL  DROPLET  GROWTH 

We  present  here  scaled  formulas  for  maximum  droplet  size  as  well  as 
location  contours  of  cloud  formation  and  evaporation.  For  each  l value 
for  distance  from  the  explosion  center,  we  may  approximate  that  portion 
of  the  M function  greater  than  unity  as  a parabola 

M(t,A)  - (M*(4)  - 1)  (1  - y2)  + 1 (4.1) 

where  M*(£)  is  the  maximum  of  M for  the  given  l (See  Figure  2);  and 


y - [(t-t2)  + (At/2)]/(At/2) 


(4.2) 


where  At(£)  = t2U)  - t^U)  and  t-j  and  t2  are  the  times  at  which  M passes 
through  unity. 

Integrating  the  scaled  droplet  growth  Equation  3.4,  subject  to  the 

universal  initial  condition  Equation  3.6,  over  the  time  such  that  S M > 

o — 

1.0  (i.e.,  between  t-j  x^*and  t2  £ t2)  we  find 

rL  - r2(v<*).*>  = 


(4/3)  At(JL) 


re  S ^ 

[ -O,  ,1 

P - <sth<‘>'So>  1 

ASth(*Jj  J 

L i - sthU) 

1/2 


(4.3) 


At («,)  is  plotted  in  Figure  2 along  with  Sth(£).  R2axU)  is  Plotted  in 
Figure  4 for  various  ambient 

By  using  the  exact  M functions  to  integrate  Equation  3.4  beyond 

o 

x2*U)  to  evaporation  where  R returns  to  zero,  we  obtain  the  evaporation 
time  t.j*U).  Figure  5 gives  contours  of  cloud  formation  (x^*(£)),  maximum 
growth  (t2*(0)>  and  evaporation  (t3*(£)),  for  various  values  of  ambient  relative 
humidity  percentages  ranging  from  70%  - 100%.  Note  that  as  the  relative  humidityap- 

proaches  100%,  the  cloud  once  formed  at  a given  location  l,  requires  an 
increasingly  longer  time  to  evaporate.  The  contours  may  be  read  at  constant 
t to  find  the  location  of  the  outside  and  inside  of  the  cloud,  as  well  as  the 


i 
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2 value  of  maximum  growth.  Figure  4 and  Equation  3.5  may  then  be  used 
to  find  maximum  droplet  radius. 

The  scaling  laws  may  be  used  to  find  the  limits  of  validity  for  the 
approach  taken  here.  We  have  noted  in  Section  2 that  for  E of  the  order 
of  10  kt  and  neglecting  water  condensation,  AX  results  in  errors  of  S of 
the  order  of  1%.  Since  droplet  radii  scale  as  E1,/6,  AX  scales  as  E1^2. 
Thus  for  E of  the  order  of  100  kt,  the  error  in  S would  be  of  order  3% 
which  is  probably  unacceptably  large  compared  to  a typical  supersaturation 
(S-l)  inside  the  cloud.  Thus,  for  accurate  results  larger  energies  re- 
quire a treatment  which  couples  cloud  formation  to  shock  hydrodynamics, 
the  simple  scaling  laws  presented  here  being  inadequate. 

The  universal  droplet  growth  equations  will  also  show  breakdown  for 
very  low  explosion  energies  such  that  the  length  scale  introduced  by  the 
ambient  droplet  sizes  becomes  important.  This  is  the  case  if  the  typical 

o 

values  of  R v (2.)  given  by  Equation  4.2  or  Figure  4 are  not  considerably 
larger  than  R2  = r£  t t where  rQ(m)  is  the  equilibrium  radius  of  a repre- 
sentative droplet  species  at  100%  relative  humidity.  Let  us  consider  a 
median  salt  grain  such  as  #3  of  Figure  3 which  has  r 


4.5  x 10”5  cm. 


This  is  well  above  the 
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For  E = 10  kt  this  corresponds  to  R ~ .5  x 10  . 

2 0 

typical  values  of  R which  in  turn  are  conservatively  of  the  order  of 
max  _l/3 

.01.  However,  since  R£  scales  like  E we  expect  our  treatment  to  break 

0 _3 

down  for  E less  than  about  10  kt. 
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FIGURE  CAPTIONS 
FIGURE  1: 

FIGURE  2: 

FIGURE  3: 

FIGURE  4: 

FIGURE  5: 


Universal  shocked  relative  humidity  multiplying  function,  M, 
versus  dynamic  time,  t,  for  air  parcels  at  various  dynamic 
lengths,  Z,  (curve  labels). 

Relative  humidity  threshold  for  cloud  formation,  S^(£), 
versus  dynamic  length  £;  M*(£)  = [S^U)]”^;  pulse  duration 
for  M > 1.0,  At(£) , versus  dynamic  length  Z. 

Droplet  growth  and  evaporation  curves  in  the  relative  humidity, 
S,  and  droplet  radius,  r,  plane.  S (r,  m),  the  equilibrium 
curves  are  shown  as  dashed  curves.  The  salt  mass,  m (grams), 
and  relative  number  n of  the  species  are  shown. 

o 

Square  of  the  universal  maximum  droplet  radius,  Rmax(J')>  versus 
dynamic  length,  Z,  for  various  ambient  relative  humidities. 

Contours  of  cloud  formation  time,  t^*,  time  of  maximum  growth 
x^*,  and  evaporation  time,  t^*,  versus  dynamic  length,  Z,  for 
various  ambient  relative  humidities. 
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Navy  Department 
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ATTN  Tom  James 

ATTN  Robert  D.  Sears  D/52-14 
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R&D  Associates 
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